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Abstract 

We consider transport diffusion in a stochastic billiard in a random 
tube which is elongated in the direction of the first coordinate (the 
tube axis). Inside the random tube, which is stationary and ergodic, 
non-interacting particles move straight with constant speed. Upon 
hitting the tube walls, they are reflected randomly, according to the 
cosine law: the density of the outgoing direction is proportional to 
the cosine of the angle between this direction and the normal vector. 
Steady state transport is studied by introducing an open tube seg- 
ment as follows: We cut out a large finite segment of the tube with 
segment boundaries perpendicular to the tube axis. Particles which 
leave this piece through the segment boundaries disappear from the 
system. Through stationary injection of particles at one boundary 
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of the segment a steady state with non-vanishing stationary parti- 
cle current is maintained. We prove (i) that in the thermodynamic 
limit of an infinite open piece the coarse-grained density profile inside 
the segment is linear, and (ii) that the transport diffusion coefficient 
obtained from the ratio of stationary current and effective boundary 
density gradient equals the diffusion coefficient of a tagged particle in 
an infinite tube. Thus we prove Fick's law and equality of transport 
diffusion and self-diffusion coefficients for quite generic rough (ran- 
dom) tubes. We also study some properties of the crossing time and 
compute the Milne extrapolation length in dependence on the shape 
of the random tube. 

Keywords: cosine law, Knudsen random walk, random medium, self- 
diffusion coefficient, transport diffusion coefficient, random walk in 
random environment 

AMS 2000 subject classifications: 60K37. Secondary: 37D50, 
60J25 

1 Introduction 

Diffusion in stationary states may be encountered either in equilibrium, where 
no macroscopic mass or energy fluxes are present in a system of many diffus- 
ing particles, or away from equilibrium, where diffusion is often driven by a 
density gradient between two open segments of the surface that encloses the 
space in which particles diffuse. In equilibrium states, one is interested in 
the self- diffusion coefficient D se if, as given by the mean-square displacement 
(MSD) of a tagged particle. This quantity, also called tracer diffusion coef- 
ficient, can be measured using e.g. neutron scattering, NMR or direct video 
imaging in the case of colloidal particles. In gradient-driven non-equilibrium 
steady states, there is a particle flux between the boundaries which is propor- 
tional to the density gradient. This factor of proportionality is the so-called 
transport or collective diffusion coefficient D trans . 

Often these two diffusion coefficients cannot be measured simultaneously 
under concrete experimental conditions and the question arises whether one 
can infer knowledge about the other diffusion coefficient, given one of them. 
Generally, in dense systems these diffusion coefficients depend in a compli- 
cated fashion on the interaction between the diffusing particles. In the case of 
diffusion in microporous media, e.g. in zeolites, however, the mean free path 
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of the particles is of the order of the pore diameter or even larger. Then dif- 
fusion is dominated by the interaction of particles with the pore walls rather 
than by direct interaction between particles. In this dilute so-called Knudsen 
regime neither D se if nor D tra ns depend on the particle density anymore, but 
are just given by the low-density limits of these two quantities. One then 
expects D se if and D trans to be equal. This assumption is a fundamental in- 
put into the interpretation of many experimental data, see e.g. [H] for an 
overview of diffusion in condensed matter systems. 

Not long ago this basic tenet has been challenged by Monte-Carlo simu- 
lation of Knudsen diffusion in pores with fractal pore walls JT6J (TTJ [18] . The 
authors of these (and further) studies concluded that self-diffusion depends 
on the surface roughness of a pore, while transport diffusion is independent 
of it. In other words, the authors of [161 [T71 [18] argue that even in the low 
density limit, where the gas particle are independent of each other and in- 
teract only with the pore walls, D se if 7^ D tra ns, with a dependence of D se if 
on the details of the pore walls that D trans does not exhibit. This coun- 
terintuitive numerical finding was quickly questioned on physical grounds 
and contradicted by further simulations [21] which give approximate equal- 
ity of the two diffusion coefficients. These controversial results gave rise to 
a prolonged debate which finally led to the consensus that indeed both dif- 
fusion coefficients should agree for the Knudsen case [21]. It has remained 
open though whether these diffusion coefficients are generally exactly equal 
or only approximately to a degree depending on the details of the specific 
setting. 

A physical argument put forward in [22] suggests general equality. To see 
this one imagines the following gedankenexperiment. Imagine one colours in 
a equilibrium setting of many non-interacting particles some of these parti- 
cles without changing their properties. At some distance from this colouring 
region the colour is removed. Then these coloured particles experience a 
density gradient just as "normal" particles in an open system with the same 
pore walls would. Since the walls are essentially the same and the properties 
of coloured and uncoloured particles are the same, the statistical properties 
of the ensemble of trajectories remain unchanged. Hence one expects any 
pore roughness to have the same effect on diffusion, irrespective of whether 
one consider transport diffusion or self-diffusion. Notice, however, that this 
microscopic argument, while intuitively appealing, is far from rigorous. First, 
the precise conditions under which the independence of the diffusion coeffi- 
cients on the pore surface is supposed to be valid, is not specified. This is 
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more than a technical issue since one may easily construct surface properties 
leading to non-diffusive behaviour (cf. [7J |2U]). Second, there is no obvious 
microscopic interpretation or unique microscopic definition of the transport 
diffusion coefficient for arbitrary surface structures. D tra ns is a genuinely 
macroscopic quantity and a proof of equality between D tra ns and D se if (which 
is naturally microscopically defined through the asymptotic long-time be- 
haviour of the MSD) requires some further work and new ideas. One needs 
to establish that on large scales the Knudsen process converges to Brownian 
motion (which then also gives D se if). Moreover, in order to compare D tra ns and 
D S eif one needs a precise macroscopic definition of D tra ns which is independent 
of microscopic properties of the system. 

The first part of this programme is carried out in [7]. There we proved the 
quenched invariance principle for the horizontal projection of the particle's 
position using the method of considering the environment viewed from the 
particle. This method is useful in a number of models related to Markov 
processes in a random environment, cf. e.g. [HI [121 US]- The aim of this 
paper is to solve the second problem of defining D trans and proving equality 
with D se if. As in [7J we consider a random tube to model pore roughness. In 
contrast to [7], we now have to consider tubes of finite extension along the 
tube contour and introduce open segments at the ends of the tube. Doing 
this rigorously then clarifies some of the salient assumptions underlying the 
equality of D trans and D self . Naturally, since we are in the dilute gas limit, 
there is no dependence on the particle density in either of the two diffusion 
constants. This obvious point has not been controversial and will not be 
stressed below. 

We note that we define D tran s through stationary transport in an open 
system since this is accessible experimentally as well as numerically in Monte 
Carlo simulation. Indeed, in the literature that gave rise to the controversy 
that we address here, this way of defining D trans is used, albeit in a non- 
rigorous fashion. Sticking to this experimentally motivated setting we shall 
give below a precise definition that can be used to prove rigorously that under 
rather generic circumstances D trans = D se i f , which means that both diffusion 
constants depend on the pore surface in the same way. As pointed out 
above, this equality is expected from independence of the particles and the 
invariance principle for the process and its time-reversed. However, we could 
not find a general result applying here, and moreover, as it turns out, the 
proof is not entirely trivial. There are some technical difficulties to overcome 
because the quenched invariance principle of Definition 12.21 below is not very 
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"strong" (there is no uniformity assumption on the speed of convergence as 
a function of the initial conditions) and the jumps of the embedded discrete- 
time billiard are not uniformly bounded. Let us mention here that it is 
generally difficult to obtain stronger results in the above sense, since the 
corrector technique, generally used in the proof of quenched central limit 
theorems for reversible Markov processes in random environment, is still not 
sufficiently well understood. 

To further illuminate the contents of our results we point out that in a 
bulk system the equality of the self-diffusion coefficient and the transport 
diffusion coefficient for the spread of equilibrium density fluctuations in an 
infinite system may be taken for granted in the case of particles that have no 
mutual interaction. Hence another way of stating the main conclusion of our 
work is the assertion that the transport diffusion coefficient as defined here in 
a stationary far- from- equilibrium setting coincides with the usual equilibrium 
transport diffusion coefficient. 

We also address finite-size effects coming from the fact that we are deal- 
ing with diffusion in a finite, open geometry. This causes deviations from 
bulk results for first-passage-time properties if a tagged particle starts its 
motion close to one boundary. In particular, we compute the permeation 
time and the Milne extrapolation length that characterizes the survival time 
of a particle injected at a boundary. 

As a final introductory remark, it is worth noting that the case of Knudsen 
gas with the cosine reflection law (which is the model considered in this 
paper) is particularly easy to analyse because the stationary state can be 
written in an explicit form, cf. Theorem 12.81 As explained below, this is 
related to the following facts: (i) there is no interaction between particles, 
(ii) for random billiard (i.e., a motion of only one particle in a closed domain) 
with the cosine reflection law the stationary measure is quite explicit, as 
shown in [6]. Similar questions are much more complicated when the explicit 
form of the stationary state is not known. This is the general situation for 
non-equilibrium steady states. We refer to e.g. the model of [2] (a chain of 
coupled oscillators) where one resorts to a bound on the entropy production. 

This paper is organized in the following way. In Section 12.11 we define 
the infinite random tube, and then introduce the process we call random 
billiard. In Section I2.2[ we then consider a gas of independent particles with 
absorption/injection in a finite piece of the random tube, and we formulate 
our results on the stationary measure for that gas and on the transport 
diffusion coefficient. In Section [231 we go on to formulate first passage time 
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results that concern exit from and crossing of the finite tube by a tagged 
particle. The remaining part of the paper is devoted to the proof of our 
results. In Section [3] we mainly use the reversibility of the process to obtain 
several technical facts used later. In Section H] we prove the result on the 
stationary measure of the Knudsen gas in the finite tube. Section [5TT1 contains 
the proofs of the results related to the transport diffusion coefficient, and in 
Section 15.21 we prove the results related to the crossing of the finite tube. 

2 General notations and main results 

Naively the transport diffusion coefficient in tube direction x may be de- 
fined through the diffusion equation for the probability density d t P(x, t) = 
d x (D(x)d x P(x,t)), where a possible x-dependence may originate from a spa- 
tial inhomogeneity of the tube. Denote by J the particle current in the 
system; assuming stationarity with a probability density P*(x) one has J = 
D(x)d x P*(x). With fixed external densities P + at x = L and P~ at x = 
one finds by integration J = D trans $ with density gradient $ = (P + — P~)/L 
and D^ ns = l/L J Q dxD~ 1 (x). By measuring the current and the bound- 
ary densities one can thus obtain the transport diffusion coefficient without 
having to determine the local quantity D(x). This result, however, implies 
knowledge of the local coarse-grained boundary densities P^ 1 to be able to 
make any comparison with D se if. In a real experimental setting as well as 
for a given microscopic model these boundary densities P^ 1 are difficult to 
obtain. In particular, there is no well-defined prescription where precisely 
on a microscopic scale these boundary quantities should be measured. We 
circumvent the problem of computing these quantities from microscopic con- 
siderations by considering the total number of particles in the tube rather 
than local properties of the boundary region of the tube. Together with 
proving a large-scale linear density profile in a stationary open random tube, 
one may then infer the macroscopic density gradient, see the definition ([3]) 
below. Thus one obtains a macroscopic definition of the transport diffusion 
coefficient which is independent of microscopic details of the model. 
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2.1 Definitions of the random tube and the random 
billiard 

In order to fix ideas in a mathematically rigorous form we first recall some 
notations from [7]. 

Let us formally define the random tube in R , d > 2. In this paper, R d_1 
will always stand for the linear subspace of R d which is perpendicular to the 
first coordinate vector e, we use the notation || • || for the Euclidean norm 
in R d or R d_1 . For k G {d - l,d} let B(x,s) = {y G R fe : \\x - y\\ < e} be 
the open e-neighborhood of x G R fc . Define S d_1 = {y G R d : |||/|| = 1} to be 
the unit sphere in R d . Let 

S/i = {w G § d_1 : h ■ w > 0} 

be the half-sphere looking in the direction h. For x G R d , sometimes it will 
be convenient to write x = (a,u), being a the first coordinate of x and 
u G then, a — x ■ e, and we write u = Ux, being U the projector 

on Fix some positive constant M, and define 

5 = {u G R* -1 : H| < M}. (1) 

Let A be an open connected domain in or M. d . We denote by dA the 
boundary of A and by A = A U cL4 the closure of A. 

The random tube is viewed as a stationary and ergodic process u> = 
(cj a , a G M), where w a is a subset of H; cf. [7] for a more detailed definition. 
We denote by P the law of this process; sometimes we will use the shorthand 
notation ( • ) for the expectation with respect to P. With a slight abuse of 
notation, we denote also by 

u = {(a,u) Gl d : u G u a } 

the random tube itself, where the billiard lives. Intuitively, u a is the "slice" 
obtained by crossing ui with the hyperplane {a} x We will assume 

that the domain u is defined in such a way that it is an open subset of W. d , 
and that it is connected. We write also u for the closure of u. In order to 
define the random billiard correctly, following [B], throughout this paper we 
suppose that P-almost surely du) is a (d — l)-dimensional surface satisfying 
the Lipschitz condition. This means that for any x G du there exist e x > 0, 
an affine isometry 3 X : M. d — > M. d , a function f x : R d_1 — > R such that 
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• f x satisfies Lipschitz condition, i.e., there exists a constant L x > such 
that \f x (z) - f x (z')\ < L x \\z - z'\\ for all z,z'] 

• 3 x x = 0, f x (0) = 0, and 

3 x (co n B{x, e x )) = {ze B(0, e x ) : *W> > . . . , z^)}. 

Roughly speaking, Lipschitz condition implies that any boundary point can 
be "touched" by a piece of a cone which lies fully inside the tube. This 
in its turn ensures that the (discrete-time) process cannot remain in a small 
neighborhood of some boundary point for very long time; in Section 2.2 of [B] 
one can find an example of a non-Lipschitz domain where the random billiard 
behaves in an unusual way. 

We keep the usual notation dx, dv, dh, . . . for the (d — l)-dimensional 
Lebesgue measure on S (usually restricted to u a for some a) or Haar measure 
on S . We write \A\ for the fc-dimensional Lebesgue measure in case A C 
]R fc , and Haar measure in case A C S . Also, we denote by u u the (d — 1)- 
dimensional Hausdorff measure on du; since the boundary is Lipschitz, one 
obtains that v u is locally finite (cf. the proof of Lemma 3.1 in [6]). 

We assume additionally that the boundary of P-a.e. u is z/^-a.e. continu- 
ously differentiable, and we denote by TZ W C du the set of boundary points 
where du is continuously differentiable. 

To avoid complications when cutting a (large) finite piece of the infinite 
random tube, we assume that there exists a constant M such that for P- 
almost all environments u we have the following: for any x,y G u> with 
\(x — y) ■ e| < 1 there exists a path connecting x,y that lies fully inside u 
and has length at most M. 

For all x £ TZ U , let us define the normal vector n^(x) £ § d_1 pointing 
inside the domain oj. 

We say that y £ u is seen from x £ Co if there exists /i £ § d_1 and to > 
such that x + th a oj for all t £ (0, to) an d x + t^h = y. Clearly, if y is seen 
from x then x is seen from y, and we write u x -H- y" when this occurs. 

Next, we construct the Knudsen random walk (KRW) (£ n , n = 0, 1, 2, . . .), 
which is a discrete time Markov process on du, cf. Section 2.2 of [6]. It is 
defined through its transition density K: for x, y £ <9cu 

^ v ((y-x)-n ul (x))((x-y)-n u] (y)) 

K(x, y) = 7 d ^' ^ ^ ^ ' ^ ^ 
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where 7<2 = ( /§ h • edh) is the normalizing constant, and !{•} stands for 
the indicator function. This means that, being P w , E w the quenched (i.e., with 
fixed u) probability and expectation, for any x G TZ U and any measurable 
B C du we have 

P.[£, t+ i eB\£ n = x] = J K(x,y) dv"{y). 

B 

We also refer to the Knudsen random walk as the random walk with cosine 
reflection law, since it is elementary to obtain from (J2J) that the density of 
the outgoing direction is proportional to the cosine of the angle between this 
direction and the normal vector. 

Remark 2.1 In fact, in the general setting of JE/, for unbounded domains, 
one has to consider the following possibility: at some moment the particle 
chooses the outgoing direction in such a way that, moving in this direction, 
it never hits the boundary of the domain again, thus going directly to the 
infinity. However, it is straightforward to see that, since u C R x E, in our 
situation P^-a.s. this cannot happen. 

It is immediate to obtain from fl2]) that K(-, •) is symmetric (that is, 
K(x,y) = K(y,x) for all x,y G du); for both the discrete- and continuous- 
time processes this leads to some nice reversibility properties, exploited in (6j 
E]. Clearly, K depends on u as well, but we usually do not indicate this in 
the notations in order to keep them simple. Also, let us denote by K n (-, ■) 
the 77,-step transition density; clearly, one obtains that K n is symmetric too 
for any n > 1. 

Now, we define the Knudsen stochastic billiard (KSB) [(Xt, Vt),t > 0), 
which is the main object of study in this paper. First, we do that for the 
process starting on the boundary du from the point x G du. Let Xq = 
£o; £i> £2, £3, • • • be the trajectory of the random walk, and define 

n 
fe=l 

Then, for t G [r n ,r n+ i), define 

Xt = in + (£n+l — £n)]T7 ^ — jT • 

llsn+1 — QnW 
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In Proposition 2.1 of [6 J it was shown that, provided that the boundary 
satisfies the Lipschitz condition, we have r n — > oo P^-a.s., and so X t is well- 
defined for all t > 0. The quantity X t stands for the position of the particle 
at time t; since it is not a Markov process by itself, we define also the cadlag 
version of the motion direction at time t: 



observe that V t G § d 1 . Recall also another notation from [6]: for x G u, 
v G § d_1 , define (with the convention inf = oo) 



so that h x (v) is the next point where the particle hits the boundary when 
starting at the location x with the direction v. Of course, we can define also 
the stochastic billiard starting from the interior of u> by specifying its initial 
position xo and initial direction v$: the particle starts at the position x and 
moves in the direction vq with unit speed until hitting the boundary at the 
point h XQ (vo); then, the previous construction is applied, being h Xo (vo) the 
starting boundary point. We denote by P* ,v the (quenched) law of KSB in 
the tube u starting from x with the initial direction v. 

Consider the rescaled projected trajectory = s _1 / 2 X st • e of KSB. 

Definition 2.2 We say that the quenched invariance principle holds for the 
Knudsen stochastic billiard in the infinite random tube if there exists a posi- 
tive constant a such that, /or IP '-almost allu, for any initial conditions (xq, Vq) 
such that h Xo (vo) G TZ^, the rescaled trajectory weakly converges 

to the Brownian motion as s — >■ oo. 

Also, for some of our results we will have to make more assumptions on 
the geometry of the random tube. Consider the following 

Condition T. 

(i) There exists a positive constant e and a continuous function (p : R i— > M d 



lim 

e4.o 



X, 



t+e 



-x, 



e 



h x (v) = x + v inf {t > : x + tv G dcu} G dcu U {oo}, 



such that 



inf \\<p(t) — x\\ > e, lim ip(t) ■ e 




oo, lim ip(t) ■ e = oo. 



x£R d \cu 
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(ii) In the case d > 3, we assume that there exist N,r% > such that 
for all x, y G TZ^ with \{x — y) • e| < 2 there exists n < N such that 
K n (x, y) > r%. 

(iii) In the case = 2, we assume that 

sup{|(x — y) • e| : x,y £ TZ^, x •£» w} < oo P-a.s. 

Remark 2.3 From iae /act iaat udxS and -almost all points of du 
belong to TZ W , it is straightforward to obtain that for Lebesguex Haar-almost 
all (x,v) G oj x we have h x (v) G 7^ (see Lemma 3.2 (i) of JfJj/). 

Remark 2.4 In the paper [1] we prove that, if the second moment of the 
projected jump length with respect to the stationary measure for the envi- 
ronment seen from the particle is finite (which is true for d > 3, but not 
always for d = 2), then under certain additional conditions ( related to Con- 
dition T of the present paper), the quenched invariance principle holds for 
the Knudsen stochastic billiard in the infinite random tube, cf. Theorem 2.2, 
Propositions 2.1 and 2.2 of FN. Let us comment more on the above Condi- 
tion T: 

• In [1], instead of the "uniform Doblin condition" (ii), we assumed a 
more explicit (although a bit more technical) Condition P, which implies 
that (ii) holds (see Lemma 3.6 of [1]). In fact, in the proof of the 
quenched invariance principle the technical condition of [1] is used only 
through the fact that it implies the uniform Doblin condition. 

• The assumption we made for d = 2 may seem to be too restrictive. 
However, is it only a bit more restrictive that the assumption that the 
random tube does not contain an infinite straight cylinder. As it was 
shown in Proposition 2.2 of [7]/, if the random tube contains an infinite 
straight cylinder, then the averaged second moment of the projected 
jump length is infinite in dimension 2, and so the (quenched) invariance 
principle cannot be valid. 

2.2 Gas of independent particles and evaluation of the 
transport diffusion coefficient 

Now, let us introduce the notations specific to this paper. Consider a positive 
number H (which is typically supposed to be large); denote by the part 



11 




Figure 1: On the definition of finite tube V' 



of the random tube uj which lies between and H : 

V% = {zELu:z-ee [0,H]}. 

Denote also 

Fff = {16 dco : x • e G (0,#)}, 
D e = {0} x u , 
D r = {H} x u H , 

so that dV% = F% U D e U D r (see Figured]). Observe that V" H can, in fact, 
consist of several separate pieces, namely, one big piece between and H, 
and possibly several small pieces near the left and the right ends (we suppose 
that H > M, so that there could not be two or more big pieces). It can be 
easily seen that those small pieces have no influence on the definition of the 
transport diffusion coefficient; for notational convention, we still allow T>^ 
to be as described above. 

Then, we consider a gas of independent particles in Vfj, described as 
follows. There is usual reflection on F£; any particle which hits U D r , 
disappears. In addition, for a given A > 0, new particles are injected in D% 
with intensity (7rf|S d_1 |) _1 A per unit surface area. Every newly injected 
particle chooses the initial direction at random according to the cosine law. 
In other words, the injection in De is given by an independent Poisson process 
in Di x § e with intensity |S d_1 | _1 A|e • u\ dxdu. 

Remark 2.5 The choice of the cosine law for the injection of new parti- 
cles is justified by Theorem 2.9 of |2J/.- for the KSB in a finite domain, the 
long-run empirical law of intersection with a (d — 1) -dimensional manifold 
is cosine. One may think of the following situation: the random tube is con- 
nected from its left side Di to a very large reservoir containing the Knudsen 
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gas in the stationary regime; then, the particles cross Di with approximately 
cosine law (at least on the time scale when the density of the particles in 
the big reservoir remains unaffected by the outflow through the tube). In 
Section^ (proof of Theorem \2.8\) we use this kind of argument to obtain a 
rigorous characterization of the steady state of this gas. 

We now consider this gas in the stationary regime. Let H[ aj &] '■= K b] x 2, 
and let Ai(a,b) be the mean number of particles in Vfj D £[ a ,6b m a fixed 
environment u. 

In Theorem 12.61 below we shall see that there exists a constant $ such 
that 



lim limsup max 



^(i - 1/2) 



H/m 



m 



which means that, after coarse-graining, the particle density profile is asymp- 
totically linear. The above quantity $ is called the (rescaled) density gradient. 

We define also the current J 1 ^ as the mean number of particles absorbed 
in D r per unit of time, and let the rescaled current be defined as 



J = lim HJ 



_ff— s-oo 



Then, consistently with the discussion in the beginning of this section, the 
transport diffusion coefficient D trans is defined by 



D t 



J 



(3) 



Now, suppose that the quenched invariance principle with constant a holds 
for the stochastic billiard. Our goal is to prove that D trans is equal to the 
self- diffusion coefficient D so if := 6~ 2 /2. To this end, we prove the following 
two results. First, we prove that the coarse-grained density profile is indeed 
linear: 

Theorem 2.6 Suppose that the quenched invariance principle holds. Then, 
for any e' > there exists m such that F-a.s. 



lim sup max 



M{^,^^) A(; -1/2) 



H/m 



ni 



l w o|> 



< e' 



(4) 
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Then, we calculate the limiting current: 

Theorem 2.7 Suppose that the quenched invariance principle holds with 
constant a, and assume also that Condition T holds. Then, we have W-a.s. 

Some remarks are in place that illustrate the significance of the above 
theorems. Theorem EL6] means that d = A(|a>o|) j and using also Theorem 12 .71 
we obtain that D trans = D se if. At the same time it becomes clear that such 
a statement can be true only asymptotically since in a finite open tube one 
has to expect finite size corrections of the mean particle number. These 
corrections may, in fact, depend strongly on the microscopic shape of the tube 
near the open boundaries. This implies that in experiments on real spatially 
inhomogeneous systems some care has to be taken as to what is measured as 
macroscopic density gradient. Notice that with Theorem 12.71 we also prove 
Fick's law for diffusive transport of matter in the random Knudsen stochastic 
billiard. Since the velocity of the particles does not change at collisions with 
the tube walls, mass transport is proportional to energy transport. In this 
interpretation Theorem 12.71 implies Fourier's law for heat conduction, see 
e.g. (21 [13] for recent work on other processes. 

For a function g G C[0, oo) and a e M, denote 

Pa (g) = mf{t>0:g(t)-g(0)=a}. (6) 

As mentioned in the introduction, in the proof of Theorems 12.61 and 12.71 we 
use the explicit form of the steady state for the Knudsen gas in the random 
tube with injection from one side. Let us formulate the following theorem: 

Theorem 2.8 

(%) For the Knudsen gas with absorption/injection in D r U L>£ (as before, 
with intensity ( f y ( i\$ d ~ 1 \)~ 1 )i per unit surface area) the unique stationary 
state is Poisson point process in T>^ x §> d ~ l with intensity A|S rf— 1 1 — 1 . 

(ii) For the gas with injection in Di only, the unique stationary distribution 
of the particle configuration is given by a Poisson point process in X?^ x 
§<2-i yjnii intensity measure 

A |grf-i|-ipK«),-^_ a ( X . e ) < pH _ a (x ■ e)] dadudh. 
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Also, in both cases, for any initial configuration the process converges to the 
stationary state described above. 

Of course, the above result is not quite unexpected. It is well known that 
independent systems have Poisson invariant distributions (with the single 
particle invariant measure for Poisson intensity), let us mention e.g. [10] 
(Section VIII. 5) and [15]. Still, we decided to include the proof of this theo- 
rem because (as far as we know), it does not directly follow from any of the 
existing results available in the literature. 

2.3 Crossing time properties 

Let us introduce some more notations for the finite random tube. We denote 
by the set of points of u , from where the particle can reach D r by a 
path which stays within V H and set Dt := {0} x u) (see Figure and 
let V H C V H be the corresponding finite tube. Since we are going to study 
now how long a tagged particle stays inside the tube and how it crosses (i.e., 
goes to the right boundary without going back to the left boundary), the 
idea is to inject it in a place from where it can actually do it. Our interest 
is then in certain first-passage properties, in particular, the total life time of 
the particle inside T> H (i.e., the time until the particle first exits T> H ) and 
the permeation time which the particle needs to first exit T> H at the end of 
the tube segment "opposite" to that where it was injected, i.e., after crossing 
the tube. 

So, suppose that one particle is injected (uniformly) at random at 
into the tube T> H (that is, the starting location has the uniform distribution 
in De, and the direction is chosen according to the cosine law), and let us 
denote by <£# the event that it crosses the tube without going back to Dg, 
i.e., <£h = { T (D r ) < r + (Di)} (here, r and r + are, respectively, entrance and 
hitting times for the discrete-time process, see (120]) and ( 121]) for the precise 
definitions). Also, define Th to be the total lifetime of the particle, i.e., if X t 
is the location of the particle at time t, then Th = min{t > : X t G DeUD r }. 

First, we calculate the asymptotic behaviour of the quenched and an- 
nealed (averaged) expectation of Th- 

Theorem 2.9 Suppose that the quenched invariance principle holds with 
constant a . We have 





P-a.s. 



(7) 



•lj I H — 



2|d> 
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Figure 2: On the definition of Di, V^, and the event <£h (a trajectory crossing 
the tube is shown) 

lim 4( E - r ") = iUilS^KM) (l^oT 1 ) • (8) 

Observe that Condition T (i) implies that |u) | is bounded away from 0, and 
so (Iwnl -1 ) < oo. At this point we remind the reader that here and in the 

V If 

next theorem the expected "times" are actually expected lengths of flight, 
related through the corresponding times through the trivial generic relation 
length = velocity x time. In our Knudsen gas we always assume unit velocity 
v — 1 so that times can be identified with the appropriate lengths. 

To elucidate the physical significance of Theorem 12.91 we observe that for 
usual Brownian motion the expected lifetime T(zo) of particle in an interval 
[0, L] is given by T(z ) = zo(L — z )/(2D), where z is the starting position 
and D is the diffusion coefficient. So, in particular, for a particle starting at 
the boundary z$ = (or at z = L) the expected life time is 0. However, in a 
microscopic model of diffusion in a finite open system, this result cannot be 
expected to be generally valid because of a positive probability that a particle 
which starts at zq = would escape through the other boundary at L. Often 
it is found empirically that the expected life time can be approximated by 

with an effective shifted coordinate z = z + Xm and effective interval length 
L = L + 2Aj\/. The empirical shift length Xm is known as Milne extrapolation 
length pE], for a recent application to diffusion in carbon nanotubes see [22J. 



16 



From the definition (jH]) one can see that the life time of a particle starting 
at the origin z = allows for the computation of the Milne extrapolation 
length through the asymptotic relation 

li m T ( Am ) = 
l^oo L 2D 

provided the diffusion coefficient D is known. 

In a physical system the Milne extrapolation length depends on molecular 
details of the gas such as type of molecule or temperature, but in a Knudsen 
gas also on the tube surface. In our model the properties of the gas are 
encoded in the unit velocity v — 1 of the particles. Observe now that the 
quantity T(Xm) corresponds to E^Th in our setting. Hence, by identifying 
H = L and using D = 6~ 2 /2, Theorem 12.91 furnishes us with the dependence 
of the Milne extrapolation length on the tube properties through 

7 d |§ a! - 1 |(|u;o|) 

Am = oil ° 2 P " a - S -' ( 10 ) 

Z\u \ 

(Am^^I^KKI^I^I- 1 )/ 2 - ( n ) 

Interestingly, Xm depends only on very few generic properties of the random 
tube. 



The next result relies on Theorem 12.71 so we need to assume a stronger 
condition on the geometry of the tube. 

Theorem 2.10 Let us suppose that the quenched invariance principle is 
valid with a, and assume that Condition T holds. For the asymptotics of 
the probability of crossing, we have 

lim H? U [C H ] = t^t (12) 

H->oo 2\0Uq\ 

lim H(P W [<£ H ]) = (\uo\- 1 ) . (13) 

H^rOO p 2 'P x 'P 

For the quenched behaviour of the conditional expectations, we have, P-a.s. 

^MT»\£h) = ± (14) 

lim -VuVhWh}) = ' | W \ (15) 

H-Hx tl t) K^O 
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Jim -E^r H m c H }) = !}\ W \ (16) 



and for the annealed ones 



^^(r H \t H )l = ^ (17) 

lim y Ew (T H I{€ H }) = ^iS^KkolUi^or 1 ) (18) 



lim ±-E„(T H I{C H }) = I 





(1^0 


l>< 

1 'P 


(l^ol 




(1^0 


l}< 

1 'P 


(Pol 



(19) 



As one sees from Theorems 12.91 and I2.1CH all our annealed results in fact 
say that one can interchange the limit as H — >■ oo with integration with 
respect to P. We still decided to include these results (even though they 
are technically not difficult) because, in models related to random environ- 
ment, it is frequent that the annealed behaviour differs substantially from 
the quenched behaviour. 

One may find it interesting to observe that, by (TT5I) and ( [TBI) 



as^oc. 

To obtain another interesting consequence of our results, let us suppose now 
that P-a.s. the random tube is such that we have \uiq\ = \uq\. Observe that, 
by Jensen's inequality, it holds that 

<i^ i> (\uo\y > i 

(and the inequality is strict if the distribution of \u \ is nondegenerate) , so 
"roughness" of the tube makes the quantities (e^Th) and (E w (7#I{£ff})) 
increase. In other words, these quantities as well as the Milne correlation 
length are minimized on the tubes with constant section (which, by the way, 
do not have to be necessarily "straight cylinders"!). 

The remaining part of the paper is devoted to the proofs of our results, 
and, as mentioned in the introduction, it is organized in the following way. In 
Section [3] we obtain several auxiliary results related to hitting of sets by the 
random billiard. In Section H] we obtain the explicit form of the stationary 
measure of the Knudsen gas in the finite tube T>^ by using the corresponding 
result from [B] about the stationary distribution of one particle in a finite 
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domain. Then, in Section 15. 1\ we apply the results of Sections E] and @] 
to obtain the explicit form of the transport diffusion coefficient. Finally, in 
Section [5T2l we use Little's theorem to prove the results related to the crossing 
time of the random tube. 

3 Some preliminary facts: hitting times and 
estimates on the crossing probabilities 

We need first to prove several auxiliary facts for random billiard in arbitrary 
finite domains. As in [6], let T> be a bounded domain with Lipschitz and a.e. 
continuously different iable boundary. We keep the notation to denote the 
law of our processes, and we still use v w to denote the (d — l)-dimensional 
Hausdorff measure on the boundary dT>. Consider a Markov chain £ on &D, 
which has a transition density K with the property K(x, y) = K(y, x) for 
all x, y G &D. Observe that the Knudsen random walk £ has the above 
property, but we need to formulate the next results in a slightly more general 
framework, since we shall need to apply them to some other processes built 
upon £. Let us introduce the notations 



for the entrance and the hitting time of B C &D. Also, for measurable B C 
du such that < v w (B) < 00 we shall write 



so that is the law for the process starting from the uniform distribution 
on B. 

Taking advantage of the reversibility of the process £, we prove the fol- 
lowing 

Lemma 3.1 Consider two arbitrary measurable sets B, F C dV such that 
Bf]F = $. 

(i) Suppose that z^(-B), z/ w (F) e (0, +00). For any F' C F , we have 



t(B) = mm{n > : £ n G B}, 
r + (B) = min{n > 1 : £ n G B} 



(20) 
(21) 




B 



P^r {F) eF'\T(F)<T+(B)] 
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1 



*>(B)P*[t(F) < r+(B)\ 

F 

1 



'(F)P£[r(B) < r+(F)] 



P y MB)<T+(F)]dv"(y) 
Vl[r{B)<T + {F)]d^{y). (22) 



F> 



(ii) Suppose that u w (B) G (0, +oo). For any B',B" C B, we have 
J K^rHB) e £V + (5) < r(F)]d^(x) 

= / K%HB) e < r(F)] di/"(x). (23) 

One immediately obtains the following consequence of Lemma 13.11 (ii) : 

Corollary 3.2 For any B,F C dV such that B n F = and v w (B) G 
(0, +oo), we /iai>e i/ie following. 

(i) For x,y G -B, Zei us define the conditional (on the event {r + (B) < 
t(F)}) transition density Kb,f{x,v) ■' 

KIZtHb) e 5" | r+(S) < t(F)} = J K B>F {x,y)dv»(y). 

B" 

Then, we have 

P*[r+(£) < T(F)]K B , F (x,y) = Pl[r + (B) < T {F)]K B , F {y,x), 

that is, the random walk conditioned to return to B without hitting F 
is reversible with the reversible measure v B F defined by 

dv^ 

-^f(x) = P^(B)<r(F)]. 

(ii) In particular (take F = in the previous part) the random walk ob- 
served at the moments of successive visits to B is reversible with the 
reversible measure . 
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Proof of Lemma \3. 11 Abbreviate for the moment U :— T> \ [B U F). First, 
write using the fact that K is symmetric 



n=l 



B {/"-i 

x / du u (x n )K(x ,Xi) . . .K(x n -i,x n ) 



F 

oo 



— > / / di/"(z„_ 1 )...di/ ,, (a;i) 

™ _i f t/™- 1 

x / dz/ w (x )i?(a; n ,a; n _i) . . .A"(a?i,ar ) 



^ P :[r(i?)=n,r + (F)>n] 



p£[t(5)<t+(F)]. 



v"(B)' 
Then, similarly 

^U F) eF'\T(F)<T+(B)] 

j oo 

= Pg[r(F)<r + (B)] g Pglr(F) = G F ' r+(5) > n] 

u-(F)pP[t(B)<t+(F)]^ i J v«{B) J K X) K 1] 

x y rfz/ w (x n )K(x ,Xi) . . . ^(rr n _i,rr n ) 

F' 

1 



u»(F)pF[r(B) < t+(F)} 
so f )22|) is proved. 



P^[r( J B)<r + (F)]^( 2/ ), 
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Let us prove (|2"3"I) . Analogously to the previous computation, we write 
J P*[£ r+(B) G 5",r+( J B) < r(F)]d^(x) 

B' 

/oo 
d^(x)J2K[^r +( B) G B",T+(B)=n,r(F) > n] 

B' n=1 
co „ „ 

= J2 I / ^(^o • • • ^(^-i) 

x J dis u (x n )K(x Q , X!) ... K(x n -i, x n ) 

B" 

CO „ „ 

= J2jd^(x n ) I ^(^...^(xi) 

x J du UJ (x Q )K(x Tl ,x n _i) . . .K(x 1 ,x ) 

B> 

= J K&HB) G B',r + (B) < r(F)}dv»(x), 

B" 

and f!23p is proved. This concludes the proof of Lemma 13.11 □ 
Next, we recall the Dirichlet's principle: 

Proposition 3.3 Consider B, F C dV with B n F = and v"{B) G 
(0, +oo), and denote h(x) = P^[t(F) < t(B)] (so that, in particular, h(x) = 
for all x G B and h(x) = 1 for all x G F). Define 

U = {h: h(x) G [0, 1], h(x) = for all x G B, h(x) = 1 for all x G F}. 

Then 

2v UJ {B)?^ ) [r{F) < t+(B)] = S{h, h) = mm£{h, h), (24) 

where 

E(h,h)= [ K{x,y){h{x)-h{y)fdv w {x)dv u {y). (25) 



(02?) 



2 
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Proof. For the proof, we refer to the discrete case, e.g. Proposition 3.8 in [T], 
and observe that the proof applies to the space-continuous case, using that, 
on general spaces, harmonicity in the analytic sense and in the probabilistic 
sense are equivalent notions by [5]. Indeed, minimizers h of the Dirichlet 
form are harmonic in the analytic sense, i.e., there are in the kernel of the 
form (see (2.10) in [5J), though the left-hand side of (24) is the value of 
£(h, h) when h is harmonic in the probabilistic sense, i.e., the expectation of 
the process at some exit time (see Theorem 2.7 in [5]) with the appropriate 
boundary conditions. □ 

Now, we go back to the Knudsen random walk in the random tube u. 
Recall that K n stands for the the n-step transition density of KRW, and 
that we have K n (x,y) = K n (y,x) for all x,y. 

Let us define for an arbitrary A C K 

F»(A) = {xedco:x-eeA}. 

In case A is an interval, say, A = [a, b), we write b) instead of F u ([a, b)). 
There is the following apriory bound on the size of the jump of the random 
billiard: there exists a constant 71 > 0, depending only on M = diam(S)/2 
and the dimension, such that for P-almost all 00 

P«[|(6 - £0) • e| > u I Co = x) < jxu-^-V, (26) 

for all x G dcu, u > 1, see formula (54) of [7j. Moreover, using (126]) . for any 
n > 1 it is straightforward to obtain that, for some 7^ > 

PMCn ~ &) • e| > u I £0 = x] < ^u-^V, (27) 

for all x G du>, u > 1 (also, without restriction of generality, we can assume 
that is nondecreasing in n). 

Now, with the help of the above formula we prove the following result: 

Lemma 3.4 For any n > 1 there exists 7^ > such that for all u > 1 and 

a GK we have 

dv"{x) J dv"{y)K n {x,y) < ^u'^. (28) 
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Proof. Abbreviate V = {a + u} x u a+u . The main idea is the following: if at 
some step the Knudsen random walk jumped from some point of oo, a + 
u) to F w [a + u, oo), it must cross V, so the probability of such a jump is the 
same as the probability of the jump to V in the semi-infinite tube with the 
boundary F w {— oo, a + u) U V. So, we obtain 



dif(x) J d^(y)K n (x,y) 

F u (- oo,a] F" [a+u,oo) 

KiCa e F"[a + u,oo))du"(x) 

F<"(-ao,a] 

n 

1 dv w (x] 



/' 6 
[J {£ fc • e > a + u, £j ■ e < a + u for all j < k} 
k=i 

F u (-oo,a] 

< f d^(x )J2 [ djf(x Q )...dif{x k .. 1 ) 

J 7 1 



k=l 

F"(-oo,a] (F^-oo.a+u))*- 1 



/ 



du UJ (x k )K(x , xt) ... K(x k -i,x k ) 

F"[a+M,oo) 

< / d^(x) I dis"(y)(K(x,y) + K 2 (x,y) + --- + K n (x,y)). 



F u (-oo,a] ^ 

By symmetry of i^, we have for any m 



F"(-oo,a) 17 V 

so Lemma [3.41 now follows from fl27l) . □ 



Let us consider a sequence of i.i.d. random variables Z\, Z 2l Z 3 , . . . with 
uniform distribution on {1,2, ...,N} (where iV is from Condition T (ii)), 
independent of everything. Also, let us define £ n := £zi+-+z n - Then, it is 
straightforward to obtain that, for any x G du and B C {y G du : — 1 < 
(y — x) ■ e < 1}, we have 

Kti e 5] > N-'n^iB) (29) 
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for some r\ > 0. Let 

N 



K(x,y) = — ^K 3 {x,y) 



N 

3=1 

be the transition density of the process (£ n ,n > 0). Observe that this process 
is still reversible with the reversible measure v w , so that K(x, y) = K(y, x) 
for all ij6 du. Similarly to [7], let us define 

Kx) = E£((£i - x) ■ e) 2 

((y-x)-e) 2 K(x,y)v"(y), (30) 



and 

b{x)=Y?M Zl -x)-ef 



((y-x)-e) 2 K(x,y)^(y). (31) 



dw 



We suppose that t(B) and r + (_B) are defined as in ([20"]) - (f2TT) but with £ 
instead of £, and let f(B) and r + (£>) be the corresponding quantities for the 
process £. 

Lemma 3.5 Suppose that B,F C <9w wif/i v u (B) e (0,oo). Moreover, as- 
sume that x ■ e < a for all x E B, and y ■ e > a + w /or all y G F fo/ 
course, the same result is valid if we assume that x ■ e < a for all x £ F, and 
y ■ e > a + u for ally E B). Then, there exist positive constants 73, 74, such 
that 

i/"(fl)P*[r(F) < r + {B)] < %u~^ + 1 J b(x) dU*{x), (32) 

F^[a,a+u] 

and 

v"(B)P*[t(F) < t+(B)} < %u~^ + ^ [ Ux) di^ix). (33) 

u\ J 

Moreover, (132]) and fl33l) are i>aM a/so m i/ie finite tube Vfj (in this case we 
assume that a > and a + u < H). 
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Proof. We keep the notation £(■,■) for the Dirichlet's form with respect to K, 
defined as in ([2"5]) . Suppose without restriction of generality that a = and 
define the function 

{0, ifx-e<0, 
1, if x • e > u, 

u^(x ■ e), if x ■ e € (0, u). 

Using Proposition 13.31 (observe that h G H) and Lemma I3.4[ we obtain 

2v»(B)?»[t(F)<t + (B)] 
< £{h,h) 

dv"(x)dv"(y)K(x,y)(h(x)-h(y)) 2 

(&j) 2 

= 2 y cf^(x) y dv"(y)K(x,y) 

F w (-oo,0) F u (u,oa) 

+ u~ 2 J dv ul (x)dv u (y)K(x,y)((y-x)-e) 2 

(F"[0,«]) 2 

+ 2u~ 2 / dv"{x) [ du u (y)K(x,y){x-e) 2 



F u [Q,u\ F w (-oo,0) 

+ 2 y di/^x) y d^(y)K(x,y){l-vT x x-ef 

F w [0,«] F u (u,oo) 

< 27 3 M"( d - 1 ) + 2w~ 2 y c/z/ w (x) y dv"{y)K{x, y){{y - x) ■ e) 2 

F u [0,u] duj 

= 27 3M ~ (d - 1) + ^ y 6(a;)di^(z), 

F w [o,a+w] 

and this proves (1321) . The proof of (1331) is completely analogous. □ 



We now work in finite tube T>%. Let us use the abbreviations U n = 
F"[n - 1, n), and K = F w [n, iJ) U £> r . Observe that, by Condition T (i), we 
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have that for some 75 G (0, +00) 

v"(U n ) > 75 (34) 

for all n and for P-a.a. u. 

To distinguish between the seconds moments of the projected jump length 
in finite and infinite tubes, we modify our notations in the following way. For 
x G dVff, let bnix) and bn(x) be the quantities defined as in fl30|) and (j3T| . 
but in the finite tube V 1 ^. Let us use the notations boo(x) and boo(x) f° r the 
corresponding quantities in the infinite tube. Now, we need an estimate on 
the integrals appearing in the right-hand sides of f l32|) and fl33|) . for the case 
of the finite tube: 

Lemma 3.6 Suppose that < S\ < S2 < 1 and assume that d > 3 and 
Condition T holds. Then, we have 

limsup— J bn{x) dif{x) < 00 F-a.s., (35) 

F^[ Sl H,s 2 H} 

and the same is valid with bu on the place ofbu- 
Proof. Let us recall some notations from [?]. Define 

& = {(co,u) : co G Q,u G duo}. 
Define the probability measure Q on & by 

dQ(u,u) = ^K^ u d^(u) dP(u), (36) 

where /Xq is the (d— 2)-dimensional Hausdorff measure on the boundary of ujq, 
kq >u is the scalar product of the normal vectors pointing inside the section and 
inside the tube (see Section 2 of [7] for details), and Z = J Q dF J„ Kq l^dfiQ (u) 
is the normalizing constant. In Lemma 3.1 of [7] it is shown that Q is the 
invariant law of the environment seen from the walker, that is 

(E u [f(0^Mn) I £0 = (0,«)]) q = (/> q . (37) 

Using also that 

■n 
k=l 
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and (1371) . it is straightforward to obtain that (&oo) Q < oo implies (&oo) Q < oo. 
So, using the notations of [7], by the ergodic theorem we obtain 

H 

J b^x) dv"{x) = J da J d^^^b^O^.v) 

F"[0,H] s 

->■ (6oo) Q as -ff — ^ oo, (38) 

a.s. and in L l , and the same with on the place of b^. Then, ( 135]) follows 
from the fact that, for all H, bn(x) < &oo(^) for all a? G i 7 ^. Now, with boo in- 
stead of &oo, the previous inequality is not necessarily valid. So, to prove ( 155]) 
for instead of 6//, consider x E du such that H^{x • e) G [si,s 2 ], and 
write (note that for all x G <9P^ we have bn{x) < i? 2 ) 



\b H (x)-b 0O (x)\ <^ 2 P^max|(a-x)-e| > ( Sl A (1 - s 2 ))H~^] 

(recall that d > 3), and then we obtain ( 1351) for bu as well. □ 
Next, we obtain a lower bound for certain escape probabilities: 

Lemma 3.7 Suppose that H/4 < n < H — 1, and m < n. Also, assume 
that d > 3 and Condition T holds. Then, there exist positive constants 77, 
7s , such that 

P^[r(K)<r + (f/ m )]>--^-, (39) 

and 



Pflr(K) < f + (D e )) > g. (40) 

Proof. Let £ be the Dirichlet form corresponding to if (cf. (1251) ). First, let 
us prove (1391 . As in Proposition I3.3[ we use the notation /i(x) = P^[f(\4) < 
observe that h(x) = for all x G C/ m and h{y) = 1 for all y & V n 
(and hence for all y G £/ n+ i)- Using this fact together with (134]) and Cauchy- 
Schwarz inequality, we write (abbreviating u := n — in) 

2if(u m )p% m [f(y n ) <f + (u m )] 

= S{h,h) 
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>J2 f du»{ Xj ) I d^(x j+1 )k(x 3 ,x j+1 )(h(x,)-h(x 3+1 )) 2 

i=Orr tt 

U+l — If f 

= ( J] ""(Urn+j)) / dv"{x ) ... / dv"{x u+1 ) 

J=0 TT TT 

U 

} j p u '(U m+j )is u '(U m+j+1 )K(xj, x j+1 )(h(xj) - h(x, J+ i)f 

j=0 

u+l , „ 

>^"V75 2 (n^(^)) _ / ^"(xo)... 
j=° ut 

/u 

j=0 

Um+u+l 

AT— 1 ~2 W+l _^ n n 

> u l^ 5 (U^(U m+j )) jdv"(x )... J dv"(x u+1 ) 



3=0 

'-' III, III — [ — t*-— J— JL 

n — m + 1 ' 

and this proves (I39I) . By denoting h(x) = P^[t(T^) < t(£^)] and writing 
2^(A)Pjlr(K)<r+(A)] 

n+1 

> 



y dv u {xj) J dv^ {xj + i)K(xj,Xj + i){h(xj) — h(xj + i)) 
+ / du CJ (x ) / rfz/ w (xi)i^(x ,Xi)(/i(xo) - h(xi)f 



D e ui 

in exactly the same way one can show (HUj) . This concludes the proof of 
Lemma 13.71 □ 

Next, we need (pointwise) estimates on the probabilities of exiting the 
tube at the left boundary: 
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Lemma 3.8 Assume Condition T andd > 3. Suppose also thatn G (~, 
and m G (0, n] . Then, there exists 79 such that for all x G U m we have 

nmt>.) < nv„)\ < M^±il . (41) 

Proof. From now on, we assume for technical reasons that m > y (in any 
case, otherwise the upper bound 1 is good enough for us). First, by Lem- 
mas 13.51 and 13.61 we obtain that 

PMf(A) < f + (f/ m )] < |. (42) 
Next, Lemma [3.71 implies that 

P^[f(V n )<f + (U m )]> n _ ( 2 +1 . (43) 

Also, from (129]) it is clear that for any x G U m we have 

K[r + (U m ) < r(D e U V n )] > P*[6 G U m ] > C 3 (44) 
for some 63 > 0. 

Now, denote <7o = f(U m ), (Tk+i = min{j > : £j G U m } to be the 
successive times when the set U m is visited. By Corollary 13.21 (i) and fj44|) . 
we obtain that, conditional on not hitting 5^11^, the process of successive 
returns to U m is reversible with the reversible density TT m (x), such that for 
all x G U m 

64 < n m (x) < C 5 

for some positive constants 64,65. Using also ( )42l) and (1131 . we obtain that 
there are constants 6*6,(77 > such that for any k 

P^[r(D e ) < f+(U m ) I r(D e U K) > cr fc ] < C " 



H 

^ m inVn) < f + (U m ) I r(D t U V n ) >a k }> ( 7 



n — m + 1 
So, we can write 

P^ m [f(A) < r(K)] = £X m [r(£<) < f(K) I f(AuK) e K-i,^]] 

k=l 
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00 C I H 

< £ nil ztt p ^ [ f u G afc 

^— ^ 07/(71 — m + 1) L 

CeC^jn-m + l) 
= 2 • (45) 

Now, the "pointwise" version of (l4"5j) is substantially more difficult to 
prove. 

Consider a sequence of i.i.d. random variables ( n G {0, 1} with 

P[Cn = 1] = ATViTs 

(recall (1291) and (1341 ). Then, one can couple the random sequences (£ n , n > 1) 
with £ = (£ n , n > 1) in such a way that when the event {( n = 1} occurs, £ n 
has the stationary distribution on ^f| n _ 1 . e ]- We denote by and E^^ the 
probability and expectation with fixed u and (, and let be the expectation 
with respect to (. One can formally define P^ in the following way. For any 
x G Ui, define the transition density R x by 



'l-N- l ri ^)R x (y) 



K(x,y), Hy<£Ui, 



Let 1Z X be the distribution on du with the density R x , and let £Yj be the uni- 
form distribution on Ui. Then, given £ n _i = 16 Ui, the law of under P^^ 
is given by 

I{Cn = l}U t + I{Cn = 0}^.. 

Also, let us define k = mm{n > 1 : ( n = 1}. 
Now, observe that 

[4- ■ e] = ■ e] on {j = R} (46) 
and, for z such that i < j, 

£ C (P:, C [l(&-&-i)-e| >«] |A = j) 
= ^(PS >f [|(| i -| i -i)-e|>u] |0 = 0) 

< ^[^] PS[l(li-^i)-e|>n] 
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< Csh-^, 
recall ( 127)) . Then, write using f )46|) and (H7)) 

^ c [max|(6-|o)-e|> S ] 



OO 



= ^ J*[ft = j] ^ (pJ* [ max | - Co) • e| > s] | « = j) 

oo 

< ^ P c [k = j]E< (P* c [there exists i < j such that 
3=1 

l(li-li-i) -e| > s/j] | k = j) 



c^-^/p^j] 

i=i 



Now, using (I48p . we have for an arbitrary x G C/ m 

pZiHDt) < f(v n )) 

= E^ ( [r(D e ) < f(V n )] 

< i^P* f [max|(x-£,) • e| < H/lQ,r(D e ) < r{V n )] 

+ P f P^Jmax|(x-40 -e| > if/ 16] 

< [max|(x-4) -e| < H/16,r{D e ) < f(V n )] + C XX H~^ . 

Let us deal with the first term in fj49|) . We have, taking advantage of 
and flUD (recall that d > 3) 

P c P^ ? [max|(x-0)-e| < H/16,f(D e ) < r(V n )] 



£>H/16 

< g ' 2( "" m+1) E^.c = 
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+ e c " (( "" m+ ; )+(m " <)) ^,fei= 

10 — 

Ci 3 (n-m + 1) 

and this concludes the proof of Lemma 13 .81 □ 

Next, we prove a result which shows that it is unlikely that a particle 
crosses the tube V H "too quickly". Suppose that one particle is injected 
(uniformly) at random at T)g into the tube T> H , and we still denote by 
the event that it crosses the tube without going back to Dg, i.e., C# = 
{r(l) r ) < r + (De)} (one can see that there is no conflict with the notation of 
Section |23|) . Also, recall that Th stands for the total lifetime of the particle 
as defined in Section |2U| i.e., if X t is the location of the particle at time t, 
then T H = min{t > : X t G D t U D r }. 

Lemma 3.9 For any e > there exists (large enough) m with the following 
property: there exists large enough Hq = Hq(u) such that for all H > Hq 

P^H^TH^m-'H 2 ]^^. (50) 

Proof. For H, m, £\ > 0, we say that x G dcu is (H, m, e^-good if 

P^[ sup \(X t -x)-e\<H/4]>l-e 1 . (51) 

Let L G Z be a large positive parameter to be specified later; for n G Z 
denote I n = F w [nL, (n + l)L); denote also 

In 1 = G /„ : x is not (H, m, £i)-good}. 

Now, consider first the case d > 3. From now on we suppose that m is 
sufficiently large to assure the following: 

P[ sup \B t \ < 1/4] >l-§, 

where B t is the standard Brownian motion and P is the corresponding prob- 
ability measure. In this case, if the invariance principle holds, then for any 
fixed E\ > every x is (H, m, £i)-good for all large enough H. Using the 



33 



monotone convergence theorem, it is straightforward to obtain that for fixed 
L,m,£i,S2 there exists large enough H such that for all H > H 

PK(/-) < e 2 \ > \. (52) 

Then, by the ergodic theorem, there exists large enough H Q such that for 
all H > H Q there exists n = n (H) such that I no C F w (iJ/4, 3F/4), and 

Now, let us consider also the event = {f(V£, no ) < r + (Z) f )} (that is, 
with respect to the process £, the particle enters Vl uo before coming back 
to Di). Then, write 

= ?^[i H ]?^[T H <m- l H 2 \l H ] 
+ pf«[£ ff ]P^fe|£ fi ]. (53) 

Now, by Lemmas 13.51 and 13. 6[ we can write for some C\ > 

max{P^[C],P^[^]}<§- (54) 

Then, from f )27|) we obtain that 

P°* [C C H | C H ] < sup P* [max ■ e - H\ > H/A] < C 2 H- {d - 1] (55) 
xeb r * v 

for some C*2 > 0. So, to complete the proof of ( 150]) . it remains to prove that 
the term pj« [7^ < m^F 2 | € H ] in fl53]) is small. 

To do this, let us recall that, by Lemma I3TTI (i), for any F' C V^ no , we 
have 

p£'[£*CW G I C H ] = {^(D^ICh])- 1 J PZ[t(D e ) < t + (V Lno )]dv«(y). 

F> 

(56) 

By Lemma [3.71 we have that for some C3 > 

(^{D^[£ H ])- 1 <C,H. (57) 
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For j > 1 denote Sj = D% U U\ U . . . U Uj . Using Lemma 13. S\ we can write 
for any y E V Lno 



Pl[r(D £ ) < T+(V Lr J} = J K{y,z)?l[r{D,) < r(V Lno )}d^(z) 

< [ K{y,z)dif{z) 



D, 



^ j 9 (Ln - j + 1) f * 



H 

_ Lno 



K{y,z)dv"{z). (58) 

So, by (1271) . in the case d > 3, we obtain from (1581) that for some positive 
constant C 4 

PllT(De)<t+(V Lno )}<^ 
and, by (!56|) . (15?]) . and the construction of n we obtain that 

G ^ I ^ < ^fe (59) 

Next, integrating (I55|) over V no \ I no , we obtain from Lemma E3] that 

_ Lno 



J P»[f (4) < r + (U ino )] A/"(y) <|£ | <fc/*(y) | 



[z)K(y,z) 



„ _, Lno 

<^^(Lno + L-i)-(^) 



if 

3=1 

< ^£ L -(d-2). 

- if 

Again using ( |56|) . ( 1571) . we obtain that 

P^[6(v ino ) G K \ /no | £*] < C 3 C 6 L-( d - 2 ). (60) 
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So, (ISTJj) and flBD]) imply that for any e 3 > there exists large enough L such 
that for all large enough H we have 

P^[^„ o) e/ no \J^|^]>l-^3. 

But then, since all x G I no \ are (if, m, £i)-good, from (l5Tj) we obtain that 

P^[Th < m- l H 2 | £ H ] < 1 - (1 - £l )(l - e 3 ). (61) 

Using ( )54|) . (|55|) . and (!6T|) in ( 153|) . we conclude the proof of (150|) in the case 
d > 3. 

Let us prove the lemma in the case d = 2. Take 

L = sup{|(x - y) ■ e| : x,y <E K u ,x & y}. 

Note that bn{x) < L 2 , so Lemma I3"31 implies that P^ £ [C//] < C-jH~ l for some 
C7 > 0. By Condition T (iii) we obtain that 

^nv Lno) eV no \l no \€ H } = 0, 

and, since for any x E Ilu 01 y £ ^[0,1/(^0 — 1)) we have K(x,y) = 0, we 
then obtain 

^r(V L n )^Ino\In 1 \^}>l-S 4 

for a small £4 > 0. The proof of f)5Up in the case d = 2 then follows in the 

same way □ 

4 On the steady state of the Knudsen gas 

In this section we prove the theorem that characterizes the stationary regime 
for the Knudsen gas in a finite tube. 

Proof of Theorem \2.8\ In order to prove item (i), we consider the process 
with absorbing/injection boundaries in both Di and D r (that is, the injection 
is given by two independent Poisson processes in x S e and D r x S(_ e ) with 
intensities |§ d_1 | _1 A|e ■ u\ dxdu in both cases). 

Fix a sequence of positive numbers Uk /* 00 such that \ut € Z for all k. 
For each k, consider a domain $^ with the following properties 

• V" H C $ fc , F£ C (A? U D r ) C 
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Figure 3: On the construction of domain 

• |$ fe | = u k ; 

• any segment ab with a G Dg U D r , b G cK&fc \ Fft has length at least 

l/(2d) 

(one may construct such a domain e.g. as shown on Figure [3]). Now, let 
us consider Xuk independent particles in By Theorem 2.4 of [6], the 
unique invariant measure of this system is product of uniform measures in 
location and direction. We are going to compare this process (observed 
only on T>fj) with the process with absorbing/injection boundaries in both 
Di and D r (naturally, we assume that the injection is with the cosine law 
and with the same intensity mentioned in Theorem 12.81 Let be the 
expectation for the above process in $^ with Xuk particles, with respect to 
the invariant measure. Also, we denote by E t the expectation with respect 
to the process with absorbing/injection boundaries in DgUD r at time t, with 
the initial configuration chosen from the Poisson point process in T) ^ x S d_1 
with intensity A | S rf 1 1 1 . 

Let ip be a function on T>^ x taking values on the interval [0, 1]. For 
a configuration rj = (xi,Vi, . . . ,x r ,v r ) in Vfj x § rf_1 (which means that we 
have r particles with positions x%, . . . , x r G and vector speeds v%, . . . , v r G 
S^ 1 ), write 

r 

ip(r)) = Ylipix^Vj). 

3=1 
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Denote also by 

ib = — / ib(x,v) dx dv 

IMIS*- 1 ! J 



the mean value of ib on T)fj x E> d 1 . 
Clearly, we have 



E m = e- x ™ j2 

= exp(Afe|(^-l)). (62) 



Also, it is straightforward to obtain that 

1 , ( i - J ¥r*'- <«> 

Since, as fc — )■ oo, the binomial distribution with parameters and |X>#|/ufc 
converges to the Poisson distribution with parameter A|X*ff|, for any ib we 
have 

lim E {k) ib{r]) = E ib{rj). (64) 

Now, let us fix t and prove that for any e > 

\E^iP(r]) - E to ib(r])\ < e (65) 

for all large enough k. For this, denote by N^(to) the total number of 
particles which entered T>fj through the right boundary D r up to time to- 
For the process with absorption/injection, an elementary calculation shows 
that N^'\t ) has Poisson distribution with parameter ('y d \'S d ~ 1 \)~ 1 \t \D r \. 
Let us suppose without restriction of generality that t < u ]!^ an d denote 

6(t ) = {(x,v) eR d x 

S(- e ) : there exists t G [0, to] suc h that x + vt G -D r }; 

observe that 0(t o ) C $ fe x §(_ e )- 

Now, a particle starting in x G \ with the direction v will cross -D r 
by time to iff (x,v) G @(to). So, it is straightforward to obtain that, for 
the process in the random variable N^(to) has the binomial distribution 
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with parameters Xuk and igi-i , which converges to the Poisson distribu- 
tion with parameter ( r y c i\Ei d ^ 1 \)^ 1 Xt \D r \ a s k — > oo. Then, conditioned on 
{N^(t ) = n}, for both processes the n entering particles to D r (seen as a 
point process on D r x §(_ e ) x [0, to]) are independent, each having density 
f(x,v,t) = £q 1 \D r \~ 1 1d\v ■ e|. Observe that the same considerations apply 
also to the particles which enter through Dg. To obtain ( 155|) . we use now 
the following coupling argument. First of all, as we already know, the initial 
configurations restricted to T>fj for both processes can be successfully coupled 
with probability that converges to 1 as k — > oo. Then, by the argument we 
just presented, the same applies for the process of particles entering through 

U D r . This shows that, with large probability, both processes can be 
successfully coupled. 

Now, combining ( |64|) with (|65|) and using the fact that a point process 
is uniquely determined by its characteristic functional (cf. e.g. Section 5.5 
of [H]), we obtain that the Poisson point process in T>°^ x S d_1 with inten- 
sity A | 1 ~ x is invariant for the Knudsen gas with absorption/injection in 
D r U D e . 

As for the convergence to the stationary state and the uniqueness, this 
follows from an easy coupling argument. Indeed, consider one process start- 
ing from the invariant measure defined above, and another process starting 
from an arbitrary (fixed) configuration. The initial particles are indepen- 
dent, but the newly injected particles are the same for both processes. Then, 
since any fixed particle will eventually disappear, the coupling time is a.s. 
finite, and so the system converges to the unique stationary state. (Using 
Theorem 2.1 of [B], with some more work one can show that, for fixed tube, 
this convergence is exponentially fast; however, we do not need this kind of 
result in the present paper.) This concludes the proof of the part (i). 

Let us prove the part (ii). Still considering the process with absorption 
and injection in D r U L>e, suppose that the particles entering through D r are 
coloured red, and the particles entering through Dg are coloured green. So, 
we need to compute the stationary measure for green particles. Using the 
(quasi) reversibility of Knudsen stochastic billiard (see Theorem 2.5 of [6]), 
we obtain that, given that there is a particle in x = (a,u) with the vector 
speed h, the probability that it is green equals 

P^ u) '- h [p- a (X-e)<p H _ a (X.e)]. 

Using also the part (i), we obtain that, for the gas with injection only in D%, 
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the stationary measure is that of Poisson point process with intensity 

X\S d - 1 \- 1 Pl a ' u) '- h [p. a {X ■ e) < p H - a (X -e)}dadudh. 

Note also that convergence and uniqueness follow from the same coupling 
argument as in part (i). This concludes the proof of Theorem 12.81 □ 

Let us observe also that Theorem 12.81 allows us to characterize the sta- 
tionary measure for Knudsen gas where the injection takes place from both 
sides, but with different intensities (which are constant on Df and D r ). We 
have 

Corollary 4.1 Consider now Knudsen gas with injection from both sides, 
with respective intensities (7d|S <i_1 |) _1 A and ^IS'* -1 !) -1 // on Df and D r 
(without restriction of generality, let us suppose that A > Then, a Poisson 
point process with intensity measure 

|S d_1 |- 1 (// + (A - fi)P { *> u) '- h [p_ a (X ■ e) < p H - a (X ■ e)]) dadudh 

is the steady state of the Knudsen gas. 

Proof. Indeed, one may imagine that particles of type 1 are injected from 
both sides with intensity (7d|§ d_1 |) _1 // and particles of type 2 are injected 
only from the left with intensity (7d|§ d_1 |) _1 (A — /i), and use Theorem 12.81 

□ " 



5 Proofs of the results on transport diffusion 
and crossing time 

5 . 1 Proof of Theorems MM and I27TI 

For integers i,j,£ > define 

R ij(g) = Hp-i(g) > pj(g)}, 
Gi, jt e(g) = l{pi(g) < £, pi(g) < p-j(g)}- 

Let be a Brownian motion with diffusion constant a, starting from the 
origin; we define (being E the expectation with respect to the probability 
measure on the space where the Brownian motion is defined) 

Rij = ERij(B^) = — - , 
40 



to be the probabilities of the corresponding events for this Brownian motion. 
Fix an integer m. For (z, h) G u x S d_1 and E\ > define 

T?(z, h) =inf {s > : lE^i^fW) - R i:J \ < e u 
\E z ' h Gi,j,m(Z ) — Gij iTn \ < Ei, 

for all i, j > such that i + j — m, and all s > So}. 

(66) 

Intuitively, T^ 1 (z, h) is the scaling factor one needs to use in order to assure 
that the rescaled (and projected on e) trajectory of the Knudsen stochastic 
billiard stays sufficiently close to the Brownian motion. 

By the portmanteau theorem, observe that, if the Knudsen stochastic 
billiard starting from (z, h) satisfies the quenched invariance principle, this 
means that for any e\ > it holds that T^(z, h) < oo. Since, for P-almost 
every oj, the invariance principle holds for a.a. starting points (z, h), we have 

J d¥ l^- 1 ]- 1 ]^]- 1 J du J dhl{T^({0,u),h) <oo} = l. 
n loo gd- 1 



By the monotone convergence theorem, we obtain that for all £i, e 2 > there 
exists t £1>£2 such that 

J dF \^ d - l \- 1 \uj Q \- 1 J du J dhI{T^((0,u),h) < t £u£2 } > 1 - e 2 . (67) 



UJ §d-l 



So, using the Ergodic Theorem, we obtain for almost all u and all H large 
enough 



W-li-l 



{(z,h) G^x S^ 1 : T^(z,h) > t E1>£2 }\ < 2e 2 H{\u \). (68) 



Then, by Theorem 12.81 we can write 

b 

M{a,b) = A|S rf 1 1 1 da J du J dhP ( ™' u) ' h [p- a {X ■ e) < p H -*(X ■ e)]. 

(69) 
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Now, let us prove that the rescaled density gradient is given by i? = 
A(|w |) . 

Proof of Theorem \2.6l Fix an arbitrary e' > and suppose that m is a (large) 
integer. Consider the quantity t m -% >m -i defined by fl67|) . and suppose that 
is large enough to assure that (recall 



i— 1 i—l 



{{z,h) eV^x^- 1 :T^(z,h) >t m - 2 , m - 2 }| <2 m - 2 H(\uj \). (70) 



Abbreviate tp := (H/m) 2 and consider any integer j e [1,^]- Sup- 
pose that ft E S d -\ z eV% are such that z ■ e e (m ~^ 1)g ], and 
2 (z, ft,) < t m -2 m -2. Then, since — > £j^ 2 -2, from fl66|) we obtain that 



P* h [p_«(X ■ e) < Ph-UX ■ e)] < PJV-Ct*-*)^) < 

= ]£ h (l-i2 m _ iJ (Z<'0)) 



^ J 1 -2 

< 



< 



m 
m 



(71) 



and 



P^b_ z . e (X ■ e) < pn-U* • e)] > P* h [p_ (m _ i+1) (#d) < ft-i^)] 

= E^(l- J R m _ J+lj _ 1 (Z^)) 



> 



m 



-2 



> 



m 
J -2 

m 



(72) 



Also, by the Ergodic Theorem, we can choose large enough so that for 
all j = 1, . . . , m 



.1 n 



m I 



r U-ljg iff] 



l^o I) 



m 

77 



\oj a \ da — (Icuo!) 

\ if 



(j-l)H 



< m -1 . (73) 



So, by (J69D, (JZOD, (TTZD, ([73D, 



m . .({m— j)H (m — j + 1)H 
H \ m m 
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^ X U X ~ m- 1 - 2m- 1 (\uj \)) 



H mm 
3 



> A— <|u; |> -Am- 1 (l + 4<|w |) )• (74) 

> 'p >' 'p 

Analogously, using (1711 instead of (1721) . we obtain 
m^f(m—j)H (m — j + l)H 



H \ m m 

< A— x 2m" 2 if ( o; ) + A— x ± (( w ) + rrT 1 

< A^-(|w |) + Am- 1 (2 + 3(|o;o|) )• (75) 

777, N 'P x 'p 

Then, we obtain (JJJ) from (174"]) and (1751) . and so the proof of Theorem 12.61 is 
concluded. □ 

At this point, let us formulate an additional result which will be used in 
Section 15.21 

Proposition 5.1 Define 

jm 
H 

M* = AIS^ 1 )- 1 J da J du J dh?^- h [p„ a {X ■ e) < p H - a (X ■ e)] 

H 

x P^> h [p H _ a (X ■ e) < p„ Q (A ■ e)], (76) 

and suppose that the quenched invariance principle holds. Then, for any 
e' > there exists m such that P-a.s. 

M > *«-l/2)( m - i+m^) <e : (77) 



lim sup max 



H/m m 
Proof. The proof is quite analogous to the proof of Theorem 12.61 □ 

Now, we calculate the limiting rescaled current. 



Proof of Theorem \2. 7[ First, we obtain an upper and a lower bounds for 
Gi.j im , where % + j = m. By e.g. the formula 1.2.0.2 of [3], we have 



t 

,2 



o 
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So, for i < m 3 / 5 

Gi,j,m > P[Pi{B^) <m\- P[p-j(B^) < Pi(B&)] 



> _ m -2/s + / ' _ J_ ds (78) 



i 2 



Also, for any i = 1, . . . , m, 



z ( % 2 



In particular, for z > m 3//5 , we obtain after some elementary computations 
that there exists a positive constant 7' such that 

1/5. 



G ij , m < 7 W/ 10 exp(-^). (80) 



Next, we employ the same strategy as in the proof of Theorem 12.61 Fix 

1 /2 _____ 

a large m, and suppose that H > mt' 2 m -% is such that ( 17D1 holds. 



Now, let y be the expected number of particles that were absorbed in 
up to time if 2 /m, in the stationary regime. Clearly, we have then = ^pj m - 



So, one can write 



H 



E^Y =A|§ d - 1 T 1 J da J du J dhP^>- h [p„ a (X ■ e) > p H -*(X • e)] 



u, 
x p(a,u),h 



# 2 



Ph- q (X ■ e) < — , p H -a(X • e) < • e) 

m 

+ EjW H ,m, (81) 

where WH, m is the mean number of particles that were injected in Dp, suc- 
cessfully crossed the tube, and then hit D r before time H 2 /m. 

Suppose that z, h are such that z ■ e 6 [ ^Ef ; ^ m ~!j' ], T™ 2 (2, /i) < 
^m- 2 ,m- 2 , T™ 2 (z,—h) < t m -2 m -2. Then, analogously to (J7TJ) and (1721) . we 
write 

P^'- ft [p_,. e (X • e) < pu-UX ■ e)] > (82) 
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V z j- h [P-UX ■ e) < Ph-UX ■ e)] < 



3 + 1 



Moreover, by (178"]) . for j < m 3//5 (recall that cp = (H/m 



m 

\2\ 



-,z,h 



H 2 



p H -z- e (X ■ e) < — , p H _ z . e (X ■ e) < p- z . e (X ■ e) 
m 

> P5 fc [Pi(#*>) < m, < p. (llH )(zW)] 

>(# 

m 

,•2 



> -m 2/5 - m 2 + 



27T(7S 3 / 2 



exp 



r 



2(X 2 s 



(is. 



(83) 



(84) 



Using (!79|) . we obtain 



H 



p H _ z . e (X ■ e) < — , p H _ z . e (X ■ e) < p- z . e (X ■ e) 



m 



— lo j— l,m— j-\-X,m\ ZJ ) 



< rrT 2 + 



J 



2nas 3 / 2 



exp 



2a 2 s 



ds. 



(85) 



Thus, using ([70]), (J73J, (EJ, OBJ), (JS5J, we obtain for some C u C 2> C a > 
(observe that, in comparison to (1741) . to estimate the product of probabilities 
in (IHTj) . we have to assume that both T™ (z, h) and T™ (z,—h) are less 
than or equal to t m -2 >m -n) 



m 

HJ% = —EajY 



H 

^77 E 



2 x§«W) -m-'-to-XW).) 



j<m 3 / 5 



2m" 2 / 5 



2tt(XS 3 / 2 



exp 



2a 2 S 



ds 



> *<M) F E ( 



. mj a/27T(js 3 / 2 



exp 



2<r 2 s 



(is— 



di-m-W-Ctm- 1 ) 
m / 
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m 

,2 



—=J- — -exp( - -l—)ds. (86) 

j<.m°/° g 

To obtain the corresponding upper bound, fix an arbitrary e > and 
suppose that m is large enough so that fl50|) of Lemma l3l3l holds for those £, m. 
The term E u WH, m of (18"T|) can be estimated in the following way: 

K,W H , m < C A —P^[€ H ,T H < m- x E 2 \ < C 4 — x |-, 
to m H 

so j^EuWh,™ < C*4£- Then, analogously to ( 1561 . using also ( 1801) . we have for 
some C 5 , C 6 > 

H J^j < A— x 4m- 2 #(U|) + C 4 e 
m \ -» J + 1 



, ill 

+ A 77 E 



if ' r to 

jr'<m 3 / 5 



7(1 



x — ((bnl) +m l )(m 2 + [ —=J- exp ( ds] 



+ A— x (to - m 3/5 ) 
H 



x— ((|a;o|) m + TO ^fm 2 + 7'm 1/10 exp f - 



2a 2 

2 



j'<m 3 / 5 o 

1/5 

+ C 4 £ + Csm- 1 + Cem 1 / 10 exp ( - . (87) 

Now, observe that 

m 

lim — / ^ — — exp ( T5~ ) 

ttwoo ^ TO y y/2icas 3 / 2 V 2(T 2 S/ 

j<m J ' 5 o 

1 

— lim \ / — exp [ — — ~ J ttz 

j<m 3 i 5 o 
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lim V — / {j/V ™ )2 exv( U/V ™ )2 )dc 



oo 1 



2 , r 2 



/cfr / rfs . — — exp ( 

o o 

1 oo 

ds f r 2 ( r 2 

— dr —j=\ — — exp 



s J V^as 1 ' 2 V 2a 2 s 



1 

ds a 2 s 

~7 x ~Y 



2 ' 



With this observation, Theorem 12.71 follows from (1551) and ( 1571) . □ 



5.2 Proof of Theorems and 12TTU1 

Observe that, since the particles are independent, the Knudsen gas in the 
finite tube Vfj can be regarded as a M/G /oo queueing system; moreover, 
using e.g. Theorem 2.1 of [6] it is straightforward to obtain that the service 
time (which is the lifetime of a newly injected particle) is a random variable 
with exponential tail. Then, let us recall the following basic identity of 
queuing theory (known as Little's theorem): 

Proposition 5.2 Suppose that A a is the arrival rate, q is the mean number 
of customers in the system, and T is the mean time a customer spends in 
the system, then T = q/A a . 

Proof. See e.g. Section 5.2 of [8J. To understand intuitively why this fact 
holds true, one may reason in the following way: by large time t, the total 
time of all the customers in the system would be (approximately) qt on one 
hand, and TA a t on the other hand. □ 

Proof of Theorem \2.9l This result almost immediately follows from Theo- 
rem [221 by using Proposition 15.21 First, for the gas of independent particles 
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the arrival rate is 

= -4!Sk (88) 

recall that the particles are injected in Di only. Then, from Theorem 12.61 it 
is straightforward to obtain that for the mean number of particles qn in the 
system, we have 

qH XH{\u \) 
hm — ■ 



-ff->oo H 2 

Then, Proposition 15.21 implies (J7J). To prove the corresponding annealed 
result, note that qn < A-ff|S| by Theorem 12.81 (ii). So, applying the bounded 
convergence theorem, we obtain (jSJ). □ 

Proof of Theorem \2.1(A First, observe that in the stationary regime the 
particles leave the system at the right boundary with rate J#, and this should 
be equal to the entrance rate A^^fC^-] of the particles which cross the tube, 
with A a from (|88|) . So, (|12|) follows from Theorem 12.71 

To prove (jl~3l) . observe that, by using Lemma [3.51 with B = and F = 
D r , we obtain that for some positive constants Ci,C 2 which do not depend 
on lo 

HP UJ [€ H ]<C 1 + ^ J b(x)du"(x). 

By (138]) . the collection of random variables (HP W [£h] , H > 1) is uniformly 
integrable, and this implies ( TT3|) . 

In order to prove (|T4T) . denote by q' H the mean number of particles in the 
stationary regime that will exit at D r . Observe that, by Theorem 12.81 (ii) 
and Proposition 15.11 



hm % = A(|w |> / x(l 



A(|w |) p 

3j ) doc — 



So, using (112"]) and Proposition 15.21 we obtain ( TI41 . The relations ( IT5l) 
and ([16]) follow from dHJ and (PJ. 

Now, observe that ( IT8|) and ( !T9|) immediately follow from f|T5|) . (|T6|) . 
and ([8]), so now it remains only to prove (jTTI) . Let o"i := r + (D^), cr^+i = 
min{m > : £ m G -D^} be the moments of successive visits to for the 
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process in the finite tube. By Corollary 13. 2 [ £ CTfc is uniformly distributed in 
for all k, and so we can write 

P^[r(£> r ) < (T k ] < kP„[£ H }. (89) 

Then, using (189 p . Lemma E2J and the fact that the random variables (Zj,j > 
1) are independent of everything, we obtain 

oo 

< P <#(A-) < r(Di), <T k -i <Z X + --- + Z f(dr) < a k \ 

k=l 

< P w [Zi H h Zj a e for all i < k and all j \ r(D r ) < a k ] 

x P w [r(/) r ) < <r k ] 

oo 
k=l 

and this implies that P w [£h-] > C±/H for some C 4 > not depending on uj. 
Since q' H < Aif|H|, one obtains (II 7p from the bounded convergence theorem. 
□ 
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